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ABSTRACT 

An onboard navigation system was developed to aid the design and evaluation of algorithms used in 
autonomous satellite navigation with Global Positioning System (GPS) data. The performance of 
the algorithms designed for a GPS Receiver/Processor Assembly (R/PA) intended for Landsat-D was 
investigated during the development phases of the GPS (four to six satellites in the constellation). 
This evaluation emphasized the effects on the orbit determination accuracy of the expected user 
clock errors, GPS satellite visibility, force model approximations, and state and covariance propaga- 
tion approximations. Results are presented giving the sensitivity of orbit determination accuracy to 
these constraints. 


*Work performed under National Aeronautics and Space Administration Contract NAS 5-24300 



INTRODUCTION 


The Navstar Global Positioning System (GPS) is a Department of Defense pro- 
gram that will provide navigation information to properly equipped users. A 
constellation of up to 24 satellites in 12-hour orbits will broadcast coded signals 
from which the user’s position can be determined. The application of GPS to 
onboard satellite navigation has been previously discussed (References 1, 2). 

As part of the evaluation of the feasibility of autonomous satellite orbit deter- 
mination using GPS, an experimental GPS Receiver/Processor Assembly (R/PA) 
will be placed on Landsat-D, and the resultant orbital solution will be compared 
to that obtained using more conventional ground-based techniques. This experi- 
ment will be conducted during the early phases of GPS, during which there will 
be four to six GPS satellites available. 

The R/PA design proposed for spacecraft applications (Reference 3) consists 
of a dual-channel receiver and a Digital Equipment Corporation (DEC) LSI-11 
processor. The R/PA measures pseudorange and delta pseudorange observa- 
tions from the GPS signals, estimates the corresponding observations using 

T 

the GPS navigation message, and uses the observation residuals in a UDU 
formulation of the extended Kalman filter (EKF) to determine the user space- 
craft’s position, velocity, clock bias and bias rate, and satellite drag coefficient. 

Simulation studies are in progress to determine the accuracy attainable with this 
use of the GPS data, to identify and evaluate the primar>^ sources of error, and 
to examine the algorithms in the proposed R/PA. As an aid to such studies, an 
onboard navigation package simulator (ONPAC) was developed on a DEC 
PDP-ll/70 computer, which has computational accuracy similar to that of the 
LSl-11, The simulator is designed for both premission planning and real-time 
analysis as well as evaluation of the GPS receiver algorithms. 

The simulator is being used in this study to determine the factors affecting the 
optimum performance of the onboard processor that are mission independent. 
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Some of these results are presented here. The topics studied include data 
editing, residual smoothing, fading of the filter memory, clock modeling, state 
process noise covariance modeling, and GPS selection. 

As an aid to the use and evaluation of GPS pseudorange and delta pseudorange 
observations, the capabilities to simulate and use these observation types were 
built into the Research and Development Goddard Trajectory Determination Sys- 
tem (R&D GTDS). GPS observation can be simulated with R&D GTDS for both 
ONPAC and R&D GTDS use. 

An overview of the steps involved in simulation and use of GPS data are sum- 
marized in Figure 1. Both truth model information and simulated data are 
passed to the ONPAC program. There, the orbit estimation is done, and the 
estimated trajectory is compared to the truth model. 

DATA SIMULATION 

The force model used in generating the true user ephemeris can be selected 
from the options available to the R&D GTDS EPHEM program (Reference 4). 
These include geopotential harmonic coefficients (up to 21-by-21), drag, solar 
radiation pressure, and perturbations from the Sun, the Moon, and the other 
planets. 

Data simulation options for parameters affecting the data accuracy are listed 
in Figure 2. It should be noted that GPS satellites can be scheduled for specific 
subsets of the total simulation time span. If a GPS satellite is scheduled for 
observations only during periods when it is not visible, it is "scheduled out" of 
the data set. The default inclination will also be a modifiable option in the 
future. 

The information passed to ONPAC is summarized in Figure 3. Because all the 
computations in ONPAC are done in the Earth-centered Eaith-fixed (ECEF) 
coordinate system, this information is in ECEF coordinates. 
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Figure 1. Overview of Analysis Approach 
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Figure 


- No EPtSOR 

- Quadratic 

- Random walk 

OBSERVATION MEASUREMENT ERRORS 
GPS CLOCK ERROR MODEL OPTIONS: 


- No ERROR 

- Constant bias, uncorrelated 

- Constant bias, correlated (i.e., all 'GPS 

CLOCKS HAVE THE SAME ERROR) 

GPS CONFIGURATION OPTIONS: 

1 TO 24 GPSs (Defaults; 5 in phase I, 12 in 
PHASE II, 24 IN PHASE III) 

- 3 ORBIT Planes 

" All orbits circular; inclination = 63 degrees; 
12-hour periods 


GPS DATA 

SPACING OPTIONS; 

- 

(from 0 TO AO) 

- At, 

(from GPS n TO GPS n-^1) 

- Atj 

(FROM last GPS IN ccnste; 


GPS SELECTION: 

- All observable 

- Geometric dilution of precision (GDOP) 

- Each GPS may be scheduled for subset(s) 

OF THE total SIMULATION TIME SPAN 

GPS EPHEMERIS ERROR OPTIONS; 

- None 

- Random constants for radial and cross-track 
(H,C); linearly increasing along-track (L) to 
A randomly selected MAXIMUM: 

— Uncorrelated 
-- Orbit-wise correlated 
-- Totally ccrrelateo 

- Sinusoidal: 

-- Input is H,C,L amplitudes, period vP), 

AND along-track RATE 1 L) 

Different ^hase offset for each GPS, 
CO.MPUTATIDN BASED CM NUMBER OF GPSS 
IN THE configuration 


2. R&D GTDS GPS Data Simulation Options 
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FOR EACH OBSERVATION; 

toBs = time of observation (including user clock offset) 

T(tl<;) = USER CLOCK OFFSET = tQgg- t|^ 

= USER CLOCK DRIFT AT 

P^Dc = "observed" OBSERVATION ( PSEUDORANGE) 

P = TRUE OBSERVATION ( PSEUDORANGE) 

= "observed" OBSERVATION (DELTA PSEUDORANGE) 

I OB S 

AfO = TRUE OBSERVATION (DELTA PSEUDORANGE) 

S^oo^ = GPS POSITION AND VELOCITY VECTORS IN ECEF 

GPo GiO pno 

COORDINATES, INCLUDING THE EFFECT OF bPS 
EPHEMERIS ERRORS 

r, V = TRUE USER POSITION AND VELOCITY VECTORS IN 

ECEF COORDINATES 

GPS SATELLITE IDENTIFICATION 

TRUTH MODEL INFORMATION PROVIDED DURING DATA GAPS 

Figure 3. Simulated Data Produced for ONPAC From R&D GTDS 
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ONPAC ESTIMATION 


T 

The ONPAC estimation is done with the UDU form of the EKF, as described 
in the mathematical specifications (Reference 5). The estimation process is 
briefly described in Figure 4. The a priori state and covariance matrix can be 
either the values from the last observation processed or the input values. The 
integration of the satellite equations of motion is done with a modified Euler 
integrator. The state transition matrix is computed with a Taylor series ap- 
proximation. Studies have demonstrated that these propagation techniques 
have sufficient accuracy for nearly circular orbits for the filter as long as the 
propagation stepsize is held small, i. e. , less than 10 seconds. 

The ONPAC state vector is given in Figure 5. The clock bias and bias rate are 
estimated in position and velocity coordinates as the clock offset and drift times 
the speed of light. The user can select all nine members as the solve-for state, 
drop the drag and estimate only eight parameters, or drop the drag and clock 
drift and estimate seven parameters. 

Parameters that can be varied in the ONPAC program are listed in Table 1. 

The force model options for the user satellite are the Earth geopotential up to 
the 5-by-5 harmonics, rotation terms, and drag. The state transition matrix 
is computed with a Taylor series approximation and has only a two-body geo- 
potential contribution plus rotation and drag terms. The effect of even further 
limitations to this force model can be studied, as can the effect of the integra- 
tor stepsize. A tunable parameter study can be done with variations of the 
process noise parameters, the fading memory smoothing factor, maximum 
values for the memory factor and the residual, and the obseiw'ation measure- 
ment noise. 

RESULTS 

As part of the autonomous orbit determination evaluation, an experimental 
R/PA will be placed on Landsat-D. The proposed Landsat-D orbit was used in 
studies of the orbit determination. 
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note: all operations are done 

IN DOUBLE PRECISION 


Figure 4. ONPAC Estimation Algorithm 
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WHERE 


X 

Y 

Z 


Cartesian position 
IN Earth-centered^ 
(ECEF) coordinates 



Clock bias 

Cartesian velocity 
IN ECEF coordinates 

Clock bias rate 
Drag coefficient 


r A 
^D '^cx 

D = 

2m 

= cross-sectional area of the 

SATELLITE 

M = MASS OF THE SATELLITE 
Cj^ = CONSTANT COEFFICIENT 


Figure 5. ONPAC State Vector 


COMPONENTS 

Earth-fixed 


COMPONENTS 
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Table 1. User Options in ON PAG 


TYPE 

OPTION 

ADJUSTABLE PARAMETERS 

Acceleration Hodel 



Geopotential 

- 

Degree & order. 2-body to 5x5 

Atmospheric drag 

On/off 

Time constant. 

EKF Algorithm 



Residual test for acceptance 

On/off 

Pmax 

Fading memory 

On/off 

/9. P2 

Process noise 

On/off 

Pa > Pb^ Pd 

SOLVE-FOR parameters 

Include/exclude 


Observation measurement noise 

drag parameter 

0-2 0-2 

State Transition Matrix 

- 

Order of approximation toaL^ 

Integrator 

- 

Step size 

1 

User clock 

- 

T iME constant. 









Initial conditions for Landsat-D and the GPS satellites are given in Figure 6. 

The GPS satellites constitute the default Phase I configuration. Since the launch 
of Landsat-D is expected during the early phases of the GPS, efforts were con- 
centrated on orbit determination using the Phase I and subsets of the Phase I 
configuration. 

Sample results are presented with four different sets of simulated data for 
October 1, 1980, 0 hours to 6 hours Universal Time (UT). The data simulation 
options used in common for these data sets are given in Figure ?• 

The ONPAC options used for four sample cases used with these data sets are 
given in Table 2. In addition, all runs were done using a 5-by-5 geopotential, 

3 

drag in the force model, a state transition matrix approximated to At , and 
a 3-second stepsize. The level of process noise used was found from tunable 
parameter studies to give the best results during periods of poor visibility. 

Figure 8 shows the root-sum- square (RSS) position error for the baseline case 
and the Phase I visibility over the 6 hours of the data span* During periods when 
the fading memory is used and four or more GPS satellites are in view, the 
RSS position error is less than 10 meters* The curve has a '’flat bottomed" 
appearance found to be characteristic of the cases when fading memory is used. 

The studies discussed here have shown that the best results occur when the 
fading memory is tuned to the periods of good GPS visibility and the process 
noise covariance to periods of poor GPS visibility. The fading memory multi- 
plies the covariance matrix and inflates the entire matrix, whereas the process 
noise is additive to certain terms of the covariance. The fading memory swamps 
the effect of the process noise when they are used together. 

Cases 2 and 3, whose RSS errors are shown in Figure 9, are done with data 
sets B and C, which include GPS clock bias errors. The GPS satellites are 
selected with the geometric dilution of precision (GDOP) procedure, which, 
when six or more GPS satellites are in view, picks for observation only those 
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LAND$AT-D INITIAL CONDITIONS 


Semimajor axis 

Eccentricity 

Inclination 

Longitude of ascending node 
Argument of perigee 
Mean anomaly 
Period 


7086.901 KILOMETERS 
0.001 

98.181 DEGREES 
35A.878 DEGREES 
180.000 DEGREES 
0.000 DEGREES 
98.956 MINUTES 


GPS CONFIGURATION AND INITIAL CONDITIONS 


PHASE I CONFIGURATION 


Inclination 
Eccentricity 
Satellites 1, 2 , 3: 

- Longitude of ascending 

- Mean anomalies 
Satellites 4. 5, 6: 

- Longitude of ascending 

- Mean anomalies 
Period 


63 DEGREES 

0.0 

NODE 120 DEGREES 

100^ 140y 180 DEGREES 

NODE 240 DEGREES 

60. 100. 140 DEGREES 
12 HOURS 


Figure 6. Initial Conditions 
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Landsat-D Force Model: 

- 8x8 GEOPOTENTIAL 

- Luni-solar perturbations 

- Drag 

- Solar radiation pressure 

Quadratic user clock error: 

- Tj = 3.3360 X 10'^ seconds 

- T2 = 3.475 X 10"^^ seconds/second 

- = 5.0 X 10"^^ seconds/second2 


Observation spacing: 

- Atj^ = 0.6 SECOND 

- At2 = 6 SECONDS 

- At^ = 6 SECONDS 

Observation standard deviations: 

- <jp = 2.0 METERS 

- = 1.7 CENTIMETERS 

Options varied: 

- GPS EPHEMERIS error MODEL 

- GPS clock bias 

- GPS CONFIGURATION AND SELECTION 


Figure 7, Data Simulation Options Used for Data 
Sets A, B, C, and D From R&D GTDS 
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Table 2. ON PAG wad Data Simulation Options Used in Sample Cases 


CASE 

FADING 

MEMORY 

PARAMETERS 

Op 

iSI 

GPS EPHEMERIS 
ERROR MODEE 

GPS CEOCK 
BIAS 

SEEECTION 

Baseline Case 

Test Case 1 
(3507) 

/5* = 0.2 

P2 = 2.0 

8.1 M 

A 

Randomly selected 

(Oj|^ CTj^) = 

(5m. 5m, 10m) 

None 

Phase L 

ALL 

Fading i%mory/GPS Clock Bias 

Test Case 2 
(AOlO) 

/3 - 0.2 

P2 = 1.05 

7.5 M 

B 

Sinusoidal 
(H.CAT) = (5m.5m. 

IOm^O.OSm/sec) 

P = 2A HOURS 

Correlated 
Tq = 3 NS 

Phase L 
GDOP 

Test Case 3 
(AOll) 

Not Used 

7.5 M 

C 

Sinusoidal 
(H.CA.L) - (5m,5m> 
10m^0.05m/sec) 

P = 12 hours 

Uncorrelated 
cr j = 3 NS 

'g 

Phase L 
GDOP 

Effect of Four in Constellation 

Test Case A 
(8102) 

iS = 0.2 
P2 - 2.0 

8.1 M 

D 

Same as C above 

None 

# 2. 3,5,6 
IN Phase I 









ROOT SUN SQUARE OF POSITION ERROR 
FOR BASELINE CASE (DATA SET A) 



PHASE I - GPS SATELLITE VISIBILITY 



Figure 8. Baseline Case and Visibility 
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ROOT SUM SQUARE OF POSITION ERROR 
FOR TEST CASE 2 (DATA SET B) 



ROOT SUM SQUARE OF POSITION ERROR 
FOR TEST CASE 3 (DATA SET C) 



Time From 0 Hours on 10/1/80 (hours) 


Figure 9, Root-Sum-Square Position Errors of Test Cases 2 and 3 
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four with the best geometric distribution. If four or fewer are visible, those 
seen are picked for observation. 

The effect of the GPS clock bias error is to increase the baseline case RSS 
position error during periods of good visibility. The fading memory option used 
in case 2 causes the RSS error to drop to the minimum value more quickly than 
in case 3 at periods of good visibility and makes the curve plotted flatter than 
that in case 3. Study of the correlated versus the uncorrelated GPS clock errors 
shows very little difference in their effects. 

Data set D was simulated using only four GPS satellites from Phase I: satel- 
lites 2 and 3 in one plane and satellites 4 and 5 in another. Figure 10 shows the 
visibility and GDOP for this data set. Test case 4, whose RSS error is shown 
in Figure 11^ was run using this data set. The RSS position error grows to 
more than 300 meters during the data gaps, and the user clock is poorly esti- 
mated when fewer than four GPS satellites are visible. Given that the error 
in the data and the GPS ephemeris is approximately 7 meters, the GDOP from 
Figure 10 would predict an error in the position determination of 35 meters or 
more when four GPS satellites are visible. The results in test case 4 are in 
the 25-through-35-meter range at times of good visibility, within the range pre- 
dicted by the GDOP. 

CONCLUSIONS 

The conclusions of the studies are given below. 

• The algorithms used are sufficient for accurate orbit determination. 

• The errors in orbit estimation are less than those predicted from 
the GDOP. 

• Accurate orbit determination is possible with only four GPS satel- 
lites in the constellation. 

• The orbit determination accuracy is limited by the GPS ephemeris 
and clock accuracies. 
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10.0 


6D0P OF GPSs # 2. 3. 5. 6 FROM 
R&D GTDS PHASE I (DATA SET D) 


Geometric 
Dilution of 
Precision 

(GDOP) 


y “ V 

5.0 I ^ ^ ^ 

0.0 1.5 3.0 A. 5 6.0 


GPS VISIBILITY FOR DATA SET D 

6.0 

A.O 

Number 

2.0 
0.0 

0.0 1.5 3.0 A. 5 6.0 

Time From 0 Hours on 10/1/80 (hours) 

Figure 10. GDOP and Visibility of Data Set D 
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• The fading memory enhances the orbit determination accuracy, 
especially when the a priori knowledge of the user clock offset is 
poor. 

All computations have been done in double precision; the effect of performing 
some operations in single precision has not yet been investigated. 

The studies discussed here have shown that the filter is not overly sensitive to 
the tunable parameters. The results presented are typical of the results 
gathered from a range of parameter values. 
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